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ABSTRACT 

We examine the accuracy of a common FIP effect diagnostic, the ratio of Ne VI to Mg VI lines in the solar transition 
region. Since the two ions have quite similar contribution functions near their maxima in equilibrium, the ratio of these 
two ions is often taken to be the abundance ratio of Ne and Mg. First we compute non-equilibrium ionization effects 
on the ratio /(Ne +5 )//(Mg +5 ) of ion fractions for a variety of simple flows through the transition region. Calculating 
the spectral line ratios for a few examples, we then show that non-equilibrium effects as well as temperature and density 
dependence must be evaluated for each line ratio used in the diagnostics. 

Subject headings: Sun: corona - Sun: abundances - Sun: solar wind 



1. Introduction 

In many regions of the solar atmosphere (including the tran- 
sition region, corona and solar wind) it appears that elements 
with a first ionization potential (FIP) greater than 10 eV are 
enhanced relative to those with FIP < 10 eV. A common diag- 
nostic to quantify this effect is the measurement of the ratios of 
emission lines of Mg VI (low FIP) and Ne VI (high FIP) (e.g. 
Jordan et al. (1998); Widing & Feldman (1989)). Both ions 
produce spectral lines in the neighborhood of 400 A (minimiz- 
ing possible calibration problems), and have very similar con- 
tribution functions in coronal equilibrium (see Fig. 1). This has 
been characterized as one of the best diagnostics for FlP-related 
abundances. 

We calculate the ion fractions for both Ne and Mg ions, e.g. 
/(Ne+ 5 ) = n(Ne +5 )/n(Ne) (Fig. 1). We then compute the ion 
ratio F = /(Ne +5 )//(Mg +5 ) in equilibrium (Fig. 2). There is 
a wide plateau near T = 4 x 10 5 K. It is this feature which is ex- 
ploited when claiming this ion ratio should be insensitive to the 
details of the structure of the corona. The density ratio depends 
on the elemental abundances of Ne and Mg. The photospheric 
abundance of Ne is poorly known, and inferred from observa- 
tions elsewhere in the solar system. However, the ratio of the 
total abundances of the two elements will alter our computed 
ratio everywhere only by a constant factor. 

It can be seen in Fig. 2 that the ratio of the two ions at 
high temperatures is very different from that at temperatures 
where the maximum in the ion fraction curves occurs, even in 
the equilibrium situation. Before interpreting ratios of spectral 
lines formed by these ions as abundance ratios it is, therefore, 
necessary to establish that the plasma is predominantly at a tem- 
perature of ~ 4 x 10 5 K. To determine the distribution of ion 
fractions, needed to interpret the line ratio, one can observe a 
series of lines from multiple Ne and Mg ions (Young & Mason 
1997). 

In the present paper we presume that the elemental abun- 
dance ratio of Ne/Mg is fixed, and examine the effect of non- 
equilibrium ionization (NEI) situations. The code used to in- 
tegrate the non-equilibrium ionization state of the gas is de- 
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Fig. 1 . — Ion fractions for Ne VI and Mg VI in coronal equilibrium. 
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Fig. 2. — Ratio F of ion fractions for coronal equilibrium, and two cooling 
flow non-equilibrium ionization scenarios (constant pressure and constant den- 
sity). 

scribed elsewhere (Edgar and Chevalier 1986; Gaetz, Edgar and 
Chevalier 1988; Esser, Edgar and Brickhouse 1998). It uses 
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tabulated rate coefficients and cooling rates from Raymond & 
Smith (1977) as updated in 1993. We compute the rates in the 
low-density limit, which may affect the dielectronic recombi- 
nation rates. The detailed shapes of the curves in Fig. 1 depend 
on what set of ionization and recombination rate coefficients 
are used, but the general conclusions presented here are inde- 
pendent of such choices. 

The relevant comparison is between the time (or distance) 
scale over which the temperature changes significantly (such as 
the temperature scale height) and the ionization or recombina- 
tion time (or distance) scale. The distance scales are given by 
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where n e is the electron density, u is the flow speed, Q on and 
a rec are ionization and recombination rate coefficients in units 
of cm 3 s _1 , respectively (Maxwellian averaged < av >, where 
o is the cross section and v is the random relative speed of elec- 
trons and ions). 
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Fig. 3. — Ratio F of ion fractions vs. fiuence (/ n e dt) for four sudden cooling 
models. Initial high temperatures T h are given in MK. 



2. Non-Equilibrium Flows 
2.1. Recombining 

As an example of a recombining flow, we look at the situa- 
tion described by Edgar and Chevalier (1986), namely a hot gas 
{Th ~ 10 6 K) permitted to cool (and recombine) due only to its 
own radiation. For T < 10 6 K, the gas cools faster than it re- 
combines, producing an over-ionized situation. This is similar 
to a microflare condition. 

We have computed two cases, and the results are plotted in 
Fig. 2. The two situations are cooling at constant pressure (iso- 
barically), and at constant density (isochorically). Note that for 
T > 2.5 x 10 5 K, while the details differ a bit, the ion ratio F is 
quite similar to that in equilibrium. However, at lower tempera- 
tures, the departure from equilibrium is quite extreme, reaching 
perhaps an order of magnitude in the isobaric case. 

The time scale for temperature change is given by the cool- 
ing time, t coo i = T /{dT /dt). Since both the recombination and 
cooling rates vary directly with density, "fluence" (J n e dt, units 
cm~ 3 s), is the natural independent variable instead of time. 
(For example, for density 10 9 cm~ 3 , a fiuence of 10 10 cm~ 3 s is 
reached in 10 s.) If the recombination time \/n e a rec becomes 
longer than the cooling time, or equivalently the recombina- 
tion coefficient becomes less than the inverse cooling fiuence, 
otrec 1 /n e t coo i, recombination will lag behind the temperature 
change. This effect can be seen in Fig. 2 for T < 2.5 x 10 s K. 

As a more extreme case of a recombining plasma, we con- 
sider a gas taken to be at coronal equilibrium at an initial high 
temperature 7], ~ 10 6 K, instantaneously reset to a cold temper- 
ature, of 10 4 K. We then follow the recombination of Ne and 
Mg through the various ionization states, as functions of the 
fiuence / n e dt. We plot F in Fig. 3, and the ionization frac- 
tions of Ne+ 5 and Mg+ 5 in Fig. 5 (downflow case). Note that in 
the recombining case, the curves for the two ions are very dif- 
ferent, which leads to a change in F by an order of magnitude 
during the recombination. This is caused by the fact that the Ne 
recombines through the +5 state sooner than the Mg does. This 
example might represent an extreme case of downflow through 
the transition region. 



It is of course true that the details of this calculation depend 
on the initial conditions, in particular on the ionization states 
of Ne and Mg in the gas. Real flows in the solar corona will 
be more complex, and one must have an idea of the history of 
the gas before it begins the downward flow through transition 
region temperatures in order to know what initial conditions are 
appropriate. However, the large variability of F in this exam- 
ple shows that downflows must be treated with caution, when 
interpreting F as an abundance ratio. (See also section 2.3.) 

2.2. Ionizing 

An extreme example of an ionizing plasma is a shock, which 
instantaneously changes the temperature from 10 4 K to of order 
10 6 K. The results are independent of the initial temperature as 
long as it is low compared to the temperature range at which 
the ions form (e.g. < 10 5 K, see Fig. 1). In the ensuing flow, we 
watch as Ne and Mg atoms ionize up beyond the +5 state, and 
compute the ion ratio. 

The ratio F is shown as a function of fiuence for four differ- 
ent shock scenarios with "post-shock" temperatures 7), of 0.5, 
1.0, 1.5, and 2.0 x 10 6 K in Fig. 4. The gas is then followed 
towards the new equilibrium situation. For 7), = 2x 10 6 K, the 
ratio F varies by about one order of magnitude, but it is rel- 
atively constant and of order unity for T/ t = 1.0 x 10 6 K. Ion 
fractions / are shown for this case in 5. 

Upflows through the transition region, modeled in the fol- 
lowing section, are similar to this situation. 

2.3. Solar Wind Flows and Downflows 

The characteristics of the solar wind flow in the critical re- 
gion, from the top of the transition region at about 10 s K to the 
corona at about 10 6 K, are not well known. As an example we 
have chosen the electron density, temperature and proton flow 
speed from a model by Hansteen, Leer, and Holzer (1997) (Fig. 
6), as this model is the only one to date that calculates the onset 
of the solar wind from the mid-chromosphere to the corona self 
consistently. In the calculations of the ion fractions shown in 
Figure 7, we have assumed that the Ne and Mg ions flow with 
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Fig. 4.— The ratio F = /(Ne+ 5 ) //(Mg+ 5 ) ion fractions vs. fluence (J n e dt) 
for four shock models. Post-shock temperatures are given in MK. 
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Density, temperature and velocity vs. height for the flow model of 
Leer, and Holzer (1997). 
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Fig. 5. — Non-equilibrium ion fractions vs. fluence, for the case of 7), = 10 6 K, 
for both ionizing (upflow) and recombining (downflow) models. 



Fig. 7. — The ratio F = /(Ne+ 5 )//(Mg+ 5 ) ion fractions for the upflow model 
of Hansteen, Leer, and Holzer (1997), and for the analogous time-reversed 
downflow model. 



the same speed as the protons. Note that this non-equilibrium 
case shows a much wider region in temperature space where F 
is of order unity. Essentially the same curve is obtained if we 
arbitrarily multiply the flow speed by 10. The equilibrium case 
(Fig. 2, solid line) is shown as a comparison. 

This scenario is somewhat similar to the shock flows pre- 
sented above (Sec 2.2), since the flow is fast enough that the 
temperature changes significantly in less than an ionization 
time: 
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This similarity is perhaps not surprising as the transition region 
is quite steep, and ends at a temperature of the order of 10 6 K, 
as in one of the shock models (dashed line in Fig. 4). 

Downflows are common in the transition region. However, 
the details such as velocity or temperature profile vs. height 
are often elusive. Here we take as an example a time-reversed 
version of the upflow model above. The ionization state of the 
gas at the beginning of the downflow is needed for the model, 
but depends on the prior history of the gas. We presume the 
gas is in coronal equilibrium at two temperatures, 0.6 and 0.8 x 



10 6 K. The resulting ion ratios F are plotted in Fig. 7. They 
start (on the right of the figure) at the equilibrium curve (the 
assumed initial condition), and move to the left toward lower 
temperatures with time. Starting from 0.8 x 10 6 K, the ion ratio 
starts at low values and comes up to nearly unity. The details are 
somewhat different from the sudden cooling model discussed 
above (Sec. 2) due to the fact that the flow decelerates towards 
lower temperatures as the density rises rapidly. This causes the 
ionization state to approach coronal equilibrium at the end of 
the flow. Nonetheless, we see that this downflow also has an 
ion ratio F significantly different from unity over the part of the 
flow where the ions are abundant. 

2.4. Emissivity Calculations 

In this section we discuss the intensity ratios of the spec- 
tral lines most commonly used for abundance ratio deter- 
mination for these two ions. The Ne VI lines are transi- 
tions from 2s2p 2 2 P 3/2 -> 2s 2 2p 2 P 3/2 ,i/2 (AA401. 93, 399.82), 
and 2s2p 22 P 1/2 ^2s 2 2p 2 P 3Al/2 (AA403.26,401.14), while 
the Mg VI lines connect 2s2p 4 4 Pi/ 2 . 3 / 2 . 5 / 2 -> 2s 2 2p 3 4 S 3/2 
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(AA399.28,400.66,403.31). Emissivity calculations were per- 
formed using the APEC code (Smith & Brickhouse 2000), 
with atomic data taken from the CHIANTI database (Landi et 
al. 1999). Energy levels for both Ne VI and Mg VI are taken 
from Martin et al. (1995). Collision and oscillator strengths 
for Ne VI are from Zhang et al. (1994) and Dankwort & Trefftz 
(1978), respectively. Those for Mg VI are taken from Bhatia & 
Mason (1980). 

The interpretation of emissivity ratios is even more compli- 
cated, as the emissivities of these lines, and their ratios, de- 
pend on both density and temperature. As examples, we plot in 
Fig. 8, the ratio Ne VI A401 . 14/Mg VI A399.28 for three of the 
above models: equilibrium ionization, isobaric cooling (with 
pressure n e T = 2.6 x 10 14 cm~ 3 K), and the Hansteen, Leer, 
and Holzer (1997) upflow model. To demonstrate the density 
dependence, the equilibrium curve is shown for two densities, 
10 7 and 10 11 cm~ 3 . Note that the region of weak temperature 
dependence (7 , w2-5x 10 5 K; e.g. Fig. 2) has more or less 
disappeared. 
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Fig. 8.— The ratio Ne VI A401.1/Mg VI A399.3 of line emissivities for coro- 
nal equilibrium at two densities, for the isobaric cooling case, and for the up- 
flow model of Hansteen, Leer, and Holzer (1997). 

These examples demonstrate that non-equilibrium effects as 
well as temperature and density dependence must be evaluated 
for each line ratio used in the diagnostic. 

3. Summary and Discussion 

The Ne+ 5 to Mg+ 5 line ratio diagnostic is generally based 
on the assumption that the ions are in ionization equilibrium 
and that the plasma is predominantly at temperatures between 
2 and 5 x 10 5 K. Below and above that range, temperature ef- 
fects are important (e.g. Fig. 2). The goal of the present paper 
was to investigate the effects of non-equilibrium situations on 
F = /(Ne +5 )/ /(Mg +5 ). If the plasma is not in equilibrium, 
assumptions regarding the plasma dynamics and prior history 
have to be made, since plasma parameters describing these flow 
dynamics are often unknown. 

To shed some light on NEI effects we have considered a 
number of idealized cases. Some of these cases, isobaric and 
isochoric cooling (Fig. 2) and ionizing shock situations (Fig. 
4) maintain an ionization balance for the two ions extremely 



close to the equilibrium values. In the presence of outflow, as 
in a solar wind situation where the flow speed in the transition 
region is small < 10 km s _1 , the ratio of the two ions actually 
remains close to 1 over a larger temperature range (Fig. 7, solid 
line). However, calculating the spectral line ratios for a few ex- 
amples shows that this effect vanishes due to the presence of a 
density dependence (Fig. 8) which exists for both equilibrium 
and NEI conditions. 

In situations resembling downflows, on the other hand, in 
which the plasma cools quickly (Figs. 3 and 7, dash and dash 
dotted lines), the deviation from equilibrium might be quite sig- 
nificant. The two cases shown in Fig. 7 demonstrate that in the 
downflowing case, the ion ratio depends very much on the up- 
per boundary condition which is not known. We have made 
the simplifying assumption that the plasma is in equilibrium 
at the temperature where it starts to downflow. This does not 
at all have to be the case. The plasma could have been out of 
equilibrium in the upflowing process, and since ionizing and re- 
combining plasmas behave differently, could depart from equi- 
librium even further. 

The calculations presented in this paper demonstrate that 
temperature and density dependence as well as non-equilibrium 
effects must be carefully evaluated for each line ratio used in the 
diagnostic. 
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